The Effect of Splayed Pins on Vortex Creep and Critical Currents 



ON 
On 
On 

Or 

o 

m 



C. J. Olson, R. T. Scalettar, G. T. Zimanyi 

Department of Physics, University of California, Davis, California 95616 

(February 1, 2008) 

We study the effects of splayed columnar pins on the vortex motion using realistic London Langevin 
simulations. At low currents vortex creep is strongly suppressed, whereas the critical current j c is 
enhanced only moderately. Splaying the pins generates an increasing energy barrier against vortex 
hopping, and leads to the forced entanglement of vortices, both of which suppress creep efficiently. 
On the other hand splaying enhances kink nucleation and introduces intersecting pins, which cut off 
the energy barriers. Thus the j c enhancement is strongly parameter sensitive. We also characterize 
the angle dependence of j c , and the effect of different splaying geometries. 

PACS numbers: 74.60.Ge,74.60.Jg 
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The enhancement of critical currents in superconduc- 
tors through irradiation by heavy ions is well established 
The ions create extended columnar defects that 
localize individual vortices much more effectively than 
naturally occurring point defects. It was suggested by 
T. Hwa et al. [|| that pinning could be further improved 
by splaying these columnar defects. In recent experi- 
ments J|-§], splayed pinning has been created directly 
through multiple-step irradiation of the sample, or indi- 
rectly by irradiation through thin foils. Up to an order 
of magnitude enhancement of the critical current relative 
to parallel columnar pins has been reported [Q. 

Columnar pins which are aligned with the field direc- 
tion allow a flux line to sit in a potential minimum in 
every layer without having to pay any elastic energy. 
There are no competing tendencies which favor point over 
columnar pinning, and hence there is a clear expectation 
that columnar pinning will greatly enhance the critical 
current, as is universally born out experimentally. 

That splay should further enhance j c seems much less 
compelling. It has been argued |3| that a vortex hopping 
between two splayed pins experiences a linearly increas- 
ing energy cost, whereas once a hop originates between 
columnar pins there is no further cost for completion 
of the jump. This leads to a suppression of motion for 
splayed pins. An additional mechanism for enhancement 
of j c that we suggest is an increase in collective pinning 
effects due to greater vortex entanglement in the pres- 
ence of splayed pins. This second mechanism may be 
more effective than the first since, although the distance 
between two splayed defects increases as one moves away 
from their point of closest approach, the separation be- 
tween the pins may be very small at this closest point. 
Thus one could argue that splaying leads to an increased 
nucleation of hopping. Furthermore, the increasing en- 
ergy cost which exists for two splayed pins can be cut off 
by the additional pins that are present. 

Set against the suggestions that splay enhances pin- 
ning is also the fact that splaying undermines the key 



feature which made columnar pinning so effective in the 
first place. Vortices must elongate in order to take advan- 
tage of the pins, though this wandering is not random, 
as for point pinning. Even accepting the increasing en- 
ergy cost for hopping and entanglement arguments, it is 
not obvious that they will dominate over the very spe- 
cial topological effectiveness of columnar pinning. The 
uncertain consequences of this set of competing effects is 
reflected experimentally in the fact that the enhancement 
of pinning by splay is sensitive to the details of the splay, 
including whether the angular distribution is Gaussian 
or bimodal, the splay angle, and the material used. 

In spite of the importance of splay for critical current 
enhancement, no realistic numerical simulations of vor- 
tices interacting with splayed pinning have been reported 
up until this time. In this paper we present the first anal- 
ysis of the effects of splayed pinning using overdamped 
molecular dynamics simulations. By examining the dis- 
tribution of lengths of pinned vortex segments between 
kinks, we explicitly verify the conjecture that the energy 
cost for vortex motion increases for intermediate length 
segments in the presence of splay, as a two pin argu- 
ment suggests. We refine this picture in a crucial way by 
demonstrating how this growth is cut off by intersections 
with additional pins. This measure at the same time 
shows an enhancement of short length segments (nucle- 
ation events). By comparing splay which is in the plane 
of the Lorentz force with splay that is orthogonal to it, we 
can also address the effect of entanglement on pinning. 
We find that the combination of entanglement and con- 
finement can lead to a suppression of creep, but only in a 
somewhat limited parameter range. Finally, we empha- 
size the need to distinguish suppression of vortex creep 
and suppression of the critical current when determining 
the efficacy of different pinning configurations. 

We conduct overdamped molecular dynamics simula- 
tions using a London-Langevin model. We refer the 
reader to |9| for details. The key feature of the approach 
is the incorporation of forms for the vortex-vortex in- 
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teraction, elastic bending, and thermal Langevin forces, 
which use experimental values for the coherence and pen- 
etration lengths, and anisotropy parameter, and hence 
allows us to simulate the system without a host of ad- 
justable parameters. In |9| we verified that our choices 
accurately reproduce the experimental phase diagram. 
The current simulations are performed at T — 77K, in- 
side the glassy phase as evidenced by the creep-like E-J 
curves below. The samples are of size 106A x 106A, con- 
taining 49 vortices extending through 80 layers. 

The pinning potential representing irradiated defects 
is modeled by short-range attractive parabolic wells of 
radius r p = 0.5A which are spatially correlated on neigh- 
boring layers to form columnar pins. The wells in dif- 
ferent layers belonging to a single column have the same 
pinning energy Up, with selected from a Gaussian 
distribution with mean U p — 0.08 and standard devia- 
tion (j p = 0.012. With this choice the vortex depinning 
transition falls near j/jo = 0.08 (where jo is the BCS de- 
pairing current), and the magnitude of the pinning force 
is on average of the same order as the elastic and Lorentz 
forces. We consider several different pinning geometries: 
parallel columnar pins, in which the columns are aligned 
with the z axis; transverse bimodal splay, in which each 
pin is tilted at angles ±9 from the z axis in the plane 
transverse to the direction of vortex motion; longitudi- 
nal bimodal splay, in which the tilt plane is spanned in- 
stead by the direction of vortex motion and the z axis; 
Gaussian splay, in which 8 is selected from a Gaussian 
distribution centered about zero and the tilt direction is 
uniformly distributed in the x—y plane; and point pins, 
which are not spatially correlated between layers. 

We first study single vortex phenomena by performing 
simulations in which the number of vortices is less than 
the number of pins, N v < N p . In Fig. @, we compare two 
samples with an equal density of columnar pins. Here 
the dimensionless resistivity, p = (E / 'j)/ Pbs is defined 
in units of pss, the Bardeen-Stephen resistivity. E is the 
electric field and j is the current density. In the first sam- 
ple the pins are aligned parallel to the z axis, and in the 
second the pins are splayed at 8 = ±5.7° from the z axis, 
transverse to the direction of vortex motion. We find a 
strong suppression of the creep of vortices by splay in the 
low current regime. While the dynamic range is limited, 
when a creep type exponential fit, E ~ exp(— l/j^), is 
performed, it is consistent with an increased value of p. 

The effect of splaying decreases with increasing cur- 
rent. The critical current j c is typically defined via a 
threshold criterion p(j c ) — p t - Choosing low thresholds 
we clearly observe an enhancement of the critical current. 
For example, using p t = 10~ 4 , j c is enhanced by 20%. At 
higher thresholds this enhancement is reduced, however, 
finally disappearing for p t > 0.05. 

Now we analyze the physical mechanisms at work in 
the presence of splay. Single vortex phenomena domi- 
nate in the dilute vortex limit, considered here. At low 



applied currents, j/jo < 0.1, vortices move between pins 
by thermally activated double kinks [|l0|-|i"3|] . For paral- 
lel columnar pins, extending an already formed double 
kink does not cost extra energy [Fig. 0(a), simulation 
image]. For splayed pins with a bimodal distribution, 
double kinks between pairs of pins tilted in opposite di- 
rections [Fig. |l|(b), simulation image] experience an in- 
creasing, or confining energy barrier, since after the nu- 
cleation of the kink the unpinned vortex segment must 
keep growing longer in the high energy region between the 
pins. On the other hand, pairs of pins tilted in the same 
direction with respect to the z axis are twice as far apart 
on average than in the columnar case, so the nucleation 
of double kinks bridging parallel pins must span twice as 
long a distance, and hence costs more energy. 

We extract the energy barrier against the spreading 
of kinks by measuring the distribution function, P(lk), 
of the lengths Ik of the vortex segments between kinks. 
If the kink energy does not depend on its length, as is 
expected for columnar pins, P(h) should be roughly uni- 
form for pins of equal depth, and exhibit slow decay with 
length for pins with differing depths. If the kink energy 
grows linearly with length, as expected for splayed pins, 
P{lk) should fall off at large Ik exponentially. 

In Fig. H we show P{lk) for j/jo — 0.065. The colum- 
nar distribution is rather uniform, whereas the splayed 
one exhibits a significant enhancement at small Ik rel- 
ative to the columnar case, followed by a rapid fall in 
the intermediate region, and a slow decay at the large Ik 
regime. The enhancement at small l^s has a simple ex- 
planation: because of the splaying, the minimal distance 
between the pins is much smaller than for the columnar 
case. Therefore the nucleation of the double kinks which 
takes place here, and the formation of short segments, is 
much enhanced. The subsequent fall of P(lk) at inter- 
mediate Ik is consistent with an exponential, supporting 
the picture of a linearly increasing potential barrier. At 
large values Ik ^ 13£, however, the decay of P(lk) is 
slowed down, and P(lk) tracks the columnar distribu- 
tion. This can be attributed to the interference of addi- 
tional pins. The high energy segment between two tilted 
pins increases only until one of the kinks reaches a third 
pin intersecting the second pin onto which the vortex is 
hopping. In the present bimodal distribution this third 
pin is parallel to the first one. Thus, as the kinks slide 
further, the length of the high energy segment remains 
unchanged, in complete analogy to the columnar case. 

To consider many-vortex phenomena, we move to the 
high density limit, N v > N p . For N v = 2N p , half of the 
vortices are directly pinned, and the rest are pinned only 
by the repulsion of their pinned neighbors. At low driving 
currents, only the latter, interstitial vortices move. The 
flowing interstitial vortices are roughly aligned with the 
applied magnetic field, whereas the pinned vortices are 
tilted. Thus the interaction between these two types of 
vortices results in their entanglement. 
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Fig. |3j shows the result of our simulations. First, the 
critical current is much reduced compared to the case 
N v < N p , due to the fact that the interstitial vortices de- 
pin much more easily than their pinned neighbors. Sec- 
ond, when j c is defined using the same threshold resis- 
tivity as before, p t = 1CP 4 , we observe a factor of 2 
enhancement in j c by splay, compared to the 20% seen 
in Fig. |l|. This strongly suggests that the forced entan- 
glement of vortices is capable of impressive additional 
enhancements of the critical current. 

Here we are in the position to check whether the 
physics of splayed pins is "adiabatically connected" to 
the columnar case. It is expected that either a single 
vortex is still sticking to a single pin, after it has been 
splayed, or that several vortices will pin to a column in 
such a way that columns are either (nearly) fully occu- 
pied, or (nearly) completely empty, forming vortex-like 
"quasi lines" ||. We tested these expectations by deter- 
mining the distribution function of percentage- wise occu- 
pation of the splayed pins. Remarkably, this distribution 
function shows sharp peaks around zero and 100% occu- 
pancy, with very suppressed values in between, verifying 
the expectations. At small angles this is mostly due to 
vortices still sticking to a single pin. 

We close this section by mentioning that the enhance- 
ment of j c occurred only for a rather limited range of 
the model's parameters; for large regions of the param- 
eter space we found either minimal effects, or the re- 
duction of j c upon splaying as the magnetic field was 
further increased. This shows the powerful influence of 
the enhanced kink nucleation and cutoff of the confining 
potential by third pins. 

Experimentally up to tenfold j c enhancements were re- 
ported |^|||, exceeding but comparable to our results. 
This enhancement is strongly dependent on the mag- 
netic field and temperature, however, often diminishing 
to small values, or even turning into a reduction instead. 
This is consistent with our finding of the importance of 
nucleation and intersecting pin effects. 

Next we return to samples with N v < N p , and explore 
the dependence of j c on the angle of splaying 9. The 
inset of Fig. || displays j c {6). For small angles 9 < 10°, 
jc increases due to the increased energy barrier to kink 
spreading with increasing angle. Around 9 ~ 10° j c ex- 
hibits a maximum. It decreases smoothly for larger an- 
gles, as the vortices cease to accommodate to the pins 
in order to stay aligned with the magnetic field |]Io| , |l3| . 
This nonmonotonic behavior of j c is consistent with ex- 
periments M. 

Experimentalists have investigated several different 
pinning configurations |^]-|| . To make contact with these 
studies, we have measured the resistivity of samples with 
equal numbers of pinning elements placed in the following 
arrangements: uncorrelated point-like pinning, Gaussian 
splayed pinning, columnar pinning, as well as longitudi- 
nal and transverse bimodal splayed pinning. 



As shown in Fig. [|, we find the highest resistivity for 
point pins, lower for columnar, and the lowest for bimodal 
transverse splay. We find that the Gaussian splay pro- 
duces an enhancement of the creep relative to columnar 
pins, which is consistent with experiments Among 
the bimodal splay configurations, transverse splay sup- 
presses creep more than longitudinal splay, again in 
agreement with experiments Q . Transverse splay is more 
effective than longitudinal, because it forces the entangle- 
ment of vortices more effectively. Also, longitudinal splay 
brings the pins closer in the direction of vortex motion, 
thus helping the nucleation of the double kinks. 

Finally we return to the current dependence. As the 
applied current increases to values approaching the de- 
pinning regime, the vortices are pinned less effectively 
and spend more time between pins. The vortex motion 
is no longer dominated by kinks, weakening the single 
vortex arguments. Also, forced entanglement is no longer 
effective at preventing vortex motion near the depinning 
regime, since when a vortex begins to move, it is increas- 
ingly likely that its forward neighbor is already moving, 
or that the push from behind is sufficient to start its 
motion. Thus instead of a forced entanglement, both 
vortices move together. For all these reasons, the creep 
suppression by splay decreases with increasing currents, 
and completely disappears in the depinning regime. 

In conclusion, we have used realistic London Langevin 
simulations to study the motion of vortices in the pres- 
ence of splayed columnar defects. We found that splay- 
ing introduces a confining potential against vortex hop- 
ping. It also enhances kink nucleation and introduces 
third pins, however, cutting off this linear potential. We 
also established the importance of forced entanglement. 
Finally we analyzed the angle dependence of the criti- 
cal current, and compared different splaying geometries. 
Several of our results compare favorably to experiments. 
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FIG. 1. Resistivity p as a function of applied current j/jo 
for columnar (filled circles) and transverse bimodal splayed 
(open squares) defects, in samples containing 49 vortices, 225 
pins, and 80 layers. The splay angle 8 — ±5.7°. A clear re- 
duction of the resistivity is seen for splayed pinning at small 
applied currents. Inset: Simulation images of individual vor- 
tices forming double kinks (a): between two columnar pins, 
(b): between two splayed pins. 




FIG. 2. Distribution of the lengths of vortex segments be- 
tween kinks, P(lk), for samples with columnar (filled circles) 
and splayed (open squares) pinning, with the same param- 
eters shown in Fig. jj] at j/jo — 0.065 (in the vortex creep 
regime). For the columnar pinning, P(lk) depends weakly on 
length, but P(lk) for splayed pinning falls off more rapidly, 
indicating that the kink energy grows with length Ik- Inset: 
Critical current j c as a function of the angle ±6 between the 
z axis and the splayed pins for the system of Fig. 1. The ini- 
tial increase of j c with 6 stops above 9 ~ 10° when the pins 
become so tilted that the vortices no longer accommodate to 
them. 
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FIG. 3. Resistivity p versus driving current j/jo for a sam- 
ple containing 49 vortices and 25 pins. Filled circles: parallel 
columnar pinning; open squares: transverse bimodal splayed 
pinning. 
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FIG. 4. Resistivity p versus driving current j/jo for sev- 
eral different pinning geometries: point pinning (plus signs), 
parallel columnar pinning (filled circles), transverse bimodal 
(open squares) and longitudinal bimodal (open right trian- 
gles) splay pinning at 6 = ±5°, and Gaussian splay pinning 
with erg = 5° (open up triangles). All correlated pinning 
produces lower resistivities than point pinning. Transverse 
bimodal splay produces a lower p than columnar pinning, but 
Gaussian splay produces a higher p. 
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